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ABSTRACT 

Using a mass-limited sample of 24/im-detected, star-forming galaxies at 0.5 < z < 1 .3, we study the mass-star 
formation rate (SFR) correlation and its tightness. The correlation is well defined (er = 0.28 dex) for disk 
galaxies (n ser sic < 1 -5), while more bulge-dominated objects often have lower specific SFRs. For disk galaxies, 
a much tighter correlation (a = 0.19 dex) is obtained if the rest-frame H-band luminosity is used instead 
of stellar mass derived from multicolor photometry. The specific SFR (sSFR) correlates strongly with rest- 
frame optical colors (hence luminosity-weighted stellar age) and also with dumpiness (which likely reflects 
the molecular gas fraction). This implies that most of the observed scatter is real, despite its low level, and 
not dominated by random measurement errors. After correcting for these differential effects a remarkably 
small dispersion remains (a — 0.14 dex), suggesting that measurement errors in mass or SFR are ^ 0.10 dex, 
excluding systematic uncertainties. Measurement errors in stellar masses, the thickening of the correlation due 
to real sSFR variations, and varying completeness with stellar mass, can spuriously bias the derived slope to 
lower values due to the finite range over which observables (mass and SFR) are available. When accounting 
for these effects, the intrinsic slope for the main sequence for disk galaxies gets closer to unity. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters — galaxies: 
structure — galaxies: high-redshift 



1. INTRODUCTION 

Star-forming galaxies obey a tight correlation between stel- 
lar mass and star formation rate (SFR), from the local Uni- 
verse (e.g., Elbaz et al. 2007; based on Brinchmann et al. 
2004; Peng et al. 2010), all the way to intermediate and at 
high redshifts (e.g., Noeske et al. 2007; Elbaz et al. 2007; 
Daddi et al. 2007; Pannella et al. 2009a; Magdis et al. 2010; 
Karim et al. 2011; Daddi et al. 2009; Stark et al. 2009; 
Lee et al. 2011; Gonzalez et al. 2010; Bouwens et al. 2011). 
With a scatter of ~ 0.3 dex at all redshifts where it has been 
measured, this relation is now known as the Main Sequence 
(MS) of star-forming galaxies. This finding has several in- 
teresting implications: (1) there is a high degree of unifor- 
mity among star-forming galaxies, and stellar mass is a cru- 
cial parameter regulating the SFR; (2) fluctuations in spe- 
cific SFRs throughout the star formation histories of actively 
star-forming galaxies are minor in all but a small number of 
outliers (e.g., Rodighiero et al. 2011); (3) at high redshifts, 
the SFR of individual galaxies must increase rapidly with 
time (Daddi et al. 2007; Renzini 2009; Peng et al. 2010; Pa- 
povich et al. 201 1). 

The existence of such a tight relation raises several ques- 
tions. First, it is often believed that stellar masses and SFRs 
for individual galaxies cannot be measured to much better 
than a factor of two precision, even in relative terms. This 
raises the question whether the small spread in the mass-SFR 
correlation is actually dominated by measurement errors, im- 
plying a potentially smaller intrinsic scatter. The question ap- 
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plies to other tight relations recently discovered for normal 
galaxies, like that between mid-IR and total IR luminosity 
(Elbaz et al. 2011, scatter ~ 0.3 dex), and between SFR and 
CO luminosity (scatter ~ 0.2 dex; Daddi et al. 2010ab; Gen- 
zel et al. 2010; Sargent et al. 2012, in preparation). Intimately 
related to this is a second issue which deals with the slope of 
the correlation. Writing SFR°<: M", it is generally found that 
a ;$ 1, but results vary from ~ 0.6 to ~ 1, depending on sam- 
ple definition, the adopted SFR indicator, and (perhaps) red- 
shift (e.g., Pannella et al. 2009a; Karim et al. 2011). Finally, 
it is generally found that the mass-SFR correlation holds for 
star-forming galaxies only, but the term "star-forming galaxy" 
is somewhat ambiguous. For example, one can refer to star- 
forming BzK samples at z ~ 2 (e.g., Daddi et al. 2007; Pan- 
nella et al. 2009a), or to blue-cloud galaxies at z = to 1 (e.g., 
Elbaz et al. 2007; Peng et al. 2010), or to Lyman break galax- 
ies at z > 3. It would be desirable to explicitly study what hap- 
pens for objectively selected and complete samples of all star- 
forming galaxies down to low levels of specific SFR (hence- 
forth, sSFR), well below the MS. To start addressing these 
questions, we will explore in this letter the role of physical 
parameters, including morphology, in driving the mass-SFR 
correlation. We assume a WMAP concordance cosmology 
and a Chabrier IMF. 

2. SAMPLE SELECTION AND MEASUREMENTS 

We study a sample of galaxies at 0.5 < z < 1.3 ex- 
tracted from the /T-band selected catalog of Daddi et al. 
(2007) in GOODS-S. We consider all star-forming galaxies 
with a S/N> 3 detection and flux > 12/iJy at 24/im from 
Spitzer+MIPS (e.g., Magnelli et al. 2011). We remove 56 
AGNs significantly detected in the 2 Ms Chandra X-ray data 
(Alexander et al. 2003). There are 599 galaxies in the sample, 
of which 70% have spectroscopic redshifts and the remain- 
der have accurate photometric redshifts from Grazian et al. 
(2006). Detailed spectral energy distributions (SEDs) are 
available based on multi-color imaging from the [/-band 
through the mid-infrared, which we use to derive rest frame 
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FIG. 1 . — The stellar mass-SFR correlation for disk galaxies (« scr sic < 1 -5, 
top) and bulge dominated galaxies (n se rsic > 1.5, bottom). Lines in the top 
panels show the best fit relation and the 1 a scatter after 3-sigma clipping of 
outliers, and are repeated in the bottom panel for comparison. Points below 
completeness are not used in the fit. 



magnitudes and colors by spline-interpolation between adja- 
cent photometry. Stellar masses are derived from SED fit- 
ting (M^ ED hereafter), using a wide range of star formation 
histories, metallicity, and allowing for dust reddening, us- 
ing the method described in Le Borgne & Rocca-Volmerange 
(2002) and as used in Elbaz et al. (2007; 2011). The Spitzer 
24/im fluxes are converted into SFRs using SED templates 
from Chary & Elbaz (2001). Analyses of Spitzer and Her- 
schel far-infrared data have shown this approach to be reli- 
able for galaxies at z < 1.3 (Magnelli et al. 2009; Elbaz et al. 
2010, 2011). The infrared-based SFR is added to the unob- 
scured component computed from the UV 1500 A luminosity 
(extrapolated from the observed photometry using the best- 
fitting SED model), without correction for dust extinction. By 
comparing the latter to the total SFR, we obtain a direct esti- 
mate of the dust attenuation at 1500A (A1500) that we use to 
correct absolute magnitudes and rest frame colors for redden- 
ing, based on the Calzetti et al. (2000) extinction law. These 
extinction-corrected magnitudes and colors are used through- 
out the following analysis. 

One aspect that we want to explore is the effect of galaxy 
morphology on the mass-SFR correlation of star-forming 
galaxies. Using GALFIT (Peng et al. 2002) we model the 
HST+ACS z-band images of each object (Giavalisco et al. 
2004, release 2.0), deriving the Sersic index (nsersic) an d half- 
light radius. It is well known that star-forming galaxies be- 
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2. — Correlation between SFR and H-band luminosity. 



come clumpier at higher redshifts (e.g., Elmegreen et al. 2004; 
Forster-Schreiberet al. 2009), hence we have also investigated 
the role of dumpiness. Galaxy dumpiness is generally de- 
fined as the fraction of light in high spatial frequency struc- 
tures (Conselice 2003; Lotz et al. 2004). Here we measure 
dumpiness (S) after subtracting the GALFIT model (G) from 
the image (/) of each galaxy: 



(1) 



The images in both the numerator and denominator are con- 
volved by the PSF, and we exclude pixels within ±1.5<7 
to reduce the noise, thus measuring residual structures only 
above this intensity threshold. The average is performed over 
the segmentation map of the galaxy, as defined by SExtrac- 
tor, which is meant to produce an intensive measurement of 
dumpiness, normalized by the spatial extent of the galaxy 
(As). We would expect that this measurement should cor- 
relate better with SFR per unit surface area rather than with 
SFR itself, and we thus also use an extensive measurement 
of dumpiness defined as S ext — S x As that should be more 
closely related to the SFR. We find that S ext is better corre- 
lated with visual dumpiness classification performed for our 
sample, dumpiness is sensitive to all sorts of structures, in- 
cluding clumps but also to the strengths of spiral arms and 
anything deviating from a single Sersic fit. 
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The mass-SFR correlation with galaxies coded according to their 
frame reddening-corrected colors (top) and dumpiness (bottom). 
Fig. 1-top. 
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FIG. 4. — The mass-SFR correlation when SFR is formally corrected for 
trends in redshift, color and dumpiness (see Table 1). Lines show the best fit 
relation and the lcr scatter after 3-sigma clipping of outliers. 



3. RESULTS 

When plotting galaxies in the mass-SFR plane (Figure [TJ, 
we find that bona-fide disk galaxies (447 objects with nsersic < 
1.5, or 72% of the parent 24/im-selected sample) display a 
clear, well-defined correlation. We adopt here a lower thresh- 
old than the canonical nsersic = 2 separation of disks/ellipticals 
in order to have a cleaner sample of disks, at the expense of 
completeness. Instead, some more bulge-dominated (nsersic > 
1.5) galaxies fall near the correlation for disk galaxies, while 
others have lower sSFRs (Figure[T}, consistent with the results 
of Wuyts et al. (2011). The latter galaxies are often visually 
classified as ellipticals and have red optical colors. We con- 
clude that the mass-SFR correlation is primarily a sequence 
of star-forming disk galaxies. Comparing to color selections, 
96% of «sersic < 1 -5 galaxies fall in the blue-cloud, while con- 
versely 73% of blue-cloud galaxies have nsersic < 1.5, imply- 
ing that a blue-clould selection of star-forming galaxies re- 
sults in a higher dispersion. Including nsersic > 1-5 galax- 
ies would also bias the slope, as completeness as a function 
of sSFR varies strongly with mass in a SFR-limited sam- 
ple (we are considering only 24/im-detected galaxies), i.e., 
a Malmquist bias. Red, bulge-dominated objects with sSFR 
more than 0.6 dex below the MS contribute ^ 6% of the to- 
tal SFR density from M* ED > 5 x 10 10 M© galaxies. This 
is likely to be an overestimate because the 24/im emission 
in some of these red spheroids, rather than from star forma- 
tion, could come from winds in evolved stars (e.g., Vega et al. 



2010) or to an AGN torus, whereas the rejection of X-ray- 
detected AGNs is strongly reducing the number of main se- 
quence galaxies at the highest stellar masses (Mullaney et al. 
2012). In the following analysis, we consider only the sample 
of 338 nsersic < 1 .5 disk galaxies satisfying redshift dependent 
mass limits which guarantee completeness down to the lower 
edge of the MS. At stellar masses higher than such mass lim- 
its, galaxies can be detected to even lower sSFRs. This can 
lead to a spurious flattening of the measured slope. 

We now explore to which extent real variations of galaxy 
properties contribute to the thickness of the mass-SFR rela- 
tion. We treat SFR as the dependent variable and present mul- 
tiple linear fits of log(SFR) as a function of log mass and other 
measured galaxy properties (P,'s): 

log(SFR) = c + Cl xPi + ...c„ x P„ (2) 

where c,'s are coefficient to be determined by the fit (see 
Table 1 for the different fits performed). Given the fairly 
wide redshift range 0.5 < z < 1.3 explored, over which the 
normalization of the mass-SFR relation changes appreciably, 
we scale all SFR measurements to a common (z) = 0.9 using 
the relation: 

log(SFR z=0 .o) = log(SFR) - c z x log^^y 

where the coefficient c z is also derived from the multilinear 
regression (Table 1). By main sequence 'slope' we mean the 
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Table 1 

Summary of multilinear fit results for log(SFR) based on Eq. 2. Each row shows a particular fit with the 
physical parameters used for that fit in column 1, the relative coefficients determined in the fit (with their s/n 
between parenthesis) in column 2, the rms between individual sfr measurements and those predicted by the 
multilinear fit in column 3, and the ope of the logarithmic sfr-m* (or sfr-l h ) relation in column 4. 



Piuamulcrs (Pj) 
(1) 


Coefficients (Cj) 
Values (S/N) 
(2) 


RMS 
dex 
(3) 


Outliers rej. 
Fraction (%) 
(4) 


Co. log Mf £Z> ,log «l+z)/1.9) 






-6.631(15.3), 0.736(1 1.4), 2.791(9.1) 


0.283 


0.89 


co.logZ.#,log (0+z)/1.9) 






-1 1 .733(29.0), i .150(3 1 .4), 0.689(3.0) 


0.186 


2.20 


c ,IogMf MP ,log ((l+z)/1.9) 






-9.1 10(18.4), 0.979(20.4), 1.416(4.8) 


0.246 


1.68 


co. log M£ MP , log ((l+z)/1.9). (£/ 


-V) re sL 




-10.147(32.4), 1.080(35.6), 0.093(0.48), -0.923(21.7) 


0.154 


2.02 


co . log Ml ED , log (( 1 +z)/ 1.9), (U - 


* / )resl 




-7.329(21.9), 0.804(24.7), 1.632(6.64), -0.969(16.2) 


0.215 


2.36 


co. log M£ MP , log ((l+z)/1.9). {U 


-V) rcsl ,logS ext 




-9.145(28.8), 0.983(32.0), 0.406(2.2), -0.956(24.4), 0.149(6.4) 


0.139 


3.03 


Const, logM£ ii:l .log((l+z)/1.9), (U - V) res , . log S cx 




-6.345(18.7).0.709(21.4), 1 .917(8.0), -0.944(16.6), 0.224(7.2) 


0.208 


1.47 



Notes: M* £D , L# and Mf M/> are in solar Units; the absolute scale of 5 ext is arbitrary 



coefficient of log(M*) (or equivalent) in the fit, whereas the 
scatter refers to the rms of the difference between the individ- 
ual measures of log(SFR) and that predicted by the multilin- 
ear fit of Eq. 2. We apply an iterative 3a clipping to remove 
strong outliers in the fit, which results in rejecting from ~ 1 
to at most ~ 3% of the galaxies. 

The mass-SFR correlation for our disk sample is found to 
have a slope of 0.74 and a scatter of 0.283 dex (Figure 1). 
This is just slightly shallower than reported by Elbaz et al. 
(2007) but steeper than that of Noeske et al. (2007) at similar 
redshifts, and comparable to the slope and scatter found by 
Rodighiero et al. (201 1) for z ~ 2 galaxies. 

The rest-frame 1 .6 /im luminosity (Lh) is sometimes used as 
a proxy for the stellar mass (e.g., Gavazzi, Pierini & Boselli 
1996; Cowie & Barger 2008), so we also tried using this in the 
fitting. Remarkably, the SFR-Lh relation is even tighter, with 
a scatter of only 0.186 dex (Table 1; Figure |2j. The AL/Lh 
ratio is affected by the light from young stars, and hence it 
is somewhat affected by the SFR. However, Ly is even more 
directly affected by star formation, and yet we find that the 
SFR-Lu relation has a scatter of 0.40 dex. When using galaxy 
colors (U — V ) res t together with Lh in the fit one obtains a rela- 
tion with even smaller scatter (0.154 dex), hence colors seem 
to be the strongest factor affecting the mass-SFR sequence. 
Binning in colors produces roughly parallel sequences in the 
mass-SFR plane (Figure |3j top). The M*/Lh ratio is also 
strongly affected by color (e.g., Pannella et al. 2009b): red- 
der galaxies have lower sSFR and higher M*/Lh, which helps 
reduce the scatter in the M* — Lh correlation compared to 
the SFR-M* relation. We thus introduce an 'empirical mass' 
M EMP by fitting the M* ED as a function of L H and (U - V ) rest 
color: 

log Mf MP = - 1 .399 + 1 .042 x log(L H ) +0.339 x (17 -V) IesU 

(3) 

similar in concept to Eq. 6-7 of Daddi et al. (2004) for BzK 
galaxies and to Eq. 1 in Bell (2008) for local galaxies (but re- 
call that here quantities are reddening corrected). The disper- 
sion between M EMP and M% ED is only 0.20 dex. We find that 
the M™ P -SFR correlation has a scatter of 0.246 dex, slightly 
smaller than the scatter in M^ ED -SFR (0.283 dex, Figure 1). A 
simple interpretation of this finding would be that, somewhat 
surprisingly, the 'empirical mass' is a more accurate stellar 
mass estimator by about 0.14 dex in relative terms (subtract- 
ing in quadrature), compared to using full SED fitting with 
a large number of galaxy population synthesis models. This 
would be possible if the large number of degrees of freedom 



(on star formation histories, etc) in the libraries of fitted tem- 
plates acts to increase the uncertainties through the various de- 
generacies of actual galaxies, perhaps because actual galaxies 
have more homogeneous star formation histories than are per- 
mitted in the SED-fitting models, leading to a smaller range 
of Mt/L ratios (at a given color) in the real world compared 
to the models, while the simple fit in Eq. 2 could be effec- 
tively equivalent to introducing priors in the SED-selection. 
Demonstrating this hypothesis would require extensive simu- 
lations that are beyond the scope of this work. We emphasize 
that Equation 2 has been empirically calibrated for our galaxy 
sample with an effective z — 0.9. It is not immediately ob- 
vious if the same relation could be applied at much lower or 
higher redshifts. 

We also find that dumpiness correlates with the SFR resid- 
uals from the average relation: clumpy galaxies tend to have 
higher sSFRs (Figure [3}bottom; Table 1), an effect which is 
significant at > 6a (Table 1) when the regression includes 5 ext 
together with reddening-corrected color. The latter fit has a 
scatter of 0.139 dex (Table 1; Figure @) and a slope of 0.98, 
which is substantially steeper than the original fit. We also 
tried using «s e rsic an d size but these parameters do not appear 
to have a measurable correlation with sSFR in our sample of 
disk galaxies. 

4. DISCUSSION 

When color and dumpiness are included as fitting parame- 
ters together with stellar mass (Figure 4; Table 1), the result- 
ing correlation ultimately reaches a scatter of about 0. 14 dex, 
compared to 0.28 for the mass-SFR correlation (Figure 1). 
This implies that most ( <L 0.24 dex) of the scatter in the orig- 
inal correlation is real, i.e., due to real variations of sSFR 
that can be traced to galaxy observables (namely, color and 
dumpiness). The combined effect of measurement errors on 
M* and SFR is thus at the level of < 0.14 dex, and likely 
smaller, as our two parameters (U — V) rest and 5 ext may not 
account for all of the physical variance. If the measurement 
errors in and SFR are similar, then each of these quan- 
tities is precise to better than 0.10 dex in relative terms, or 
about 25% in linear scale, which is quite amazing. It could 
be argued, however, that what we are comparing here is not 
a direct measurement of these quantities. For example, our 
SFR is derived from the total infrared luminosity L^ \, which 
in turn is inferred from 12^m rest-frame luminosity L\i (the 
UV contribution is generally negligible), hence our conclu- 
sions strictly apply to such an observable, and similar consid- 
erations could held forM*. In principle, the exact conversion 
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between h\i and the bolometric luminosity (Lbol) could affect 
the slope that we infer for the mass-SFR relation, although 
our adopted value is perfectly consistent with all observables 
at z ~ 1 (Elbaz et al. 201 1). It could also be possible that this 
L\2 to Lbol conversion actually reduces the observed scatter, 
so that a smaller rms is found for the tracer (e.g. L\2) with 
respect to what it is designed to trace (e.g., SFR or Lbol)- This 
might be the case if the correlation between IR8 {=L%^ m / 'Lbol) 
and sSFR (Elbaz et al. 2011) were to hold also at 12/im rest 
frame and inside the main sequence. On the other hand, for 
moderately star forming galaxies mid-IR might be better cor- 
related to SFR than Lbol (see, e.g., Calzetti et al. 2007) 

It is worth commenting on the reason for the dependence 
of sSFR on color and dumpiness. For the latter the interpre- 
tation appears to be quite straightforward given that a higher 
gas fraction is expected to be the reason for higher dumpiness 
(Bournaud, Elmegreen & Elmegreen 2007; Ceverino, Dekel 
& Bournaud 2010), and the higher gas fraction should directly 
imply higher SFRs (all other things being equal, hence higher 
sSFR). Higher sSFR obviously results in bluer colors, but age, 
metallicity, or extinction effects may also be at play. We ex- 
plored if color differences could be connected to a metallicity 
dependence of the mass-SFR relation (Mannucci et al. 2010; 
Lara-Lopez et al. 2010). However, the effect of a metallicity- 
SFR relation should depend strongly on stellar mass, so it 
seems unlikely that it could produce the series of almost par- 
allel sequences as a function of color seen in Figure [3] In 
addition, metallicity has a weak impact on optical colors of 
actively star-forming galaxies. Using synthetic galaxy spec- 
tra we find that higher metallicity actually can imply slightly 
bluer colors for young galaxies. We exclude the possibility 
that colors reflect, even in part, dust reddening, because the 
luminosities and colors used in this analysis have all been cor- 
rected for dust extinction (see Section 2). We also find that 
using A 1500 as an extra parameter would further reduce the 
scatter to 0.125 dex. However, this may be due to the fact that 
A1500 is derived in part from the same infrared data used to 
determine the SFR. 

We have demonstrated that a sizable fraction of the dis- 
persion of the SFR-M* relation for main sequence galaxies 
in our sample can be traced to intrinsic differences in color 
and dumpiness at fixed mass, as SFR correlates with these 
quantities. While the existence of a main sequence of star- 
forming galaxies indicates that star formation proceeds in a 
quasi-steady fashion in most galaxies, it is quite natural to ex- 
pect that the sSFR is also subject to up and down fluctuations, 



possibly responding to fluctuations in the gas accretion rate. 
On the other hand, the star formation process is intrinsically 
stocastic, as is the formation of clumps with enhanced SFR. 
In addition, some galaxies may sustain systematically lower 
or higher SFRs, over times comparable to the Hubble time, 
which would produce quite diverging mass growth histories 
over cosmological timescales (Renzini 2009). Addressing this 
question in detail is beyond the scope of this letter. 

Finally, this work provides hints about the possible intrin- 
sic slope of the mass-SFR correlation. Indeed, the slope on 
the stellar mass term rises from 0.74 to ~ 1 when including 
color and dumpiness or using the empirical mass in the fit 
(see Table 1). When the analysis is limited to sSFR ranges 
with high completeness, the average color of galaxies in our 
sample does not change with stellar mass, which disfavors the 
slope steepening being just the effect of spurious removal of 
mass-SFR trends. Instead, using color at least partially re- 
moves the bias due to the higher dynamic range in sSFR at 
higher masses, allowing to see galaxies with lower sSFR and 
redder colors. Also, fitting for color and dumpiness removes 
most of the (real) thickness of the mass-SFR relation, which 
can also artificially reduce the slope, as we are spanning only 
a relatively small range in both stellar mass and SFR, i.e., 
about 1.5 dex in this and in similar studies. The slope will 
also be biased by additional errors in the stellar mass, which 
is used here as the independent variable but is also affected by 
measurement errors. Therefore, our work suggests that the in- 
trinsic slope of the mass-SFR correlation is probably closer to 
unity than is suggested by simple fits neglecting these biases. 

We conclude by emphasizing that the H-band luminosity 
and the (reddening corrected) rest-frame (U — V ) res t color can 
be used to predict the 24/im emission (a measure of dust and 
PAH luminosity) with better than a 40% accuracy, which is an 
impressive demonstration of the deep connection between the 
amount of stars already present (the stellar mass) and its time 
derivative (the star formation rate). 
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